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We consider an electronic spin in a quantum dot, coupled to the surrounding nuclear spins via 
inhomogeneous antiferromagnetic hyperfine interactions and subject to a uniform field, which is 
', described by Gaudin's central spin model. We study spectral properties, the two-point correla- 

■ tion functions, and the magnetization profile in the ground state and in low-lying excited states, 

' which characterizes the structure of the cloud of nuclear spins screening the electron spin. A close 

, connection to the pair occupation probability in the BCS-model is established. Using the exact 

Bethe Ansatz solution of that model and arguments of integrability, we can distinguish between 
contributions from purely classical physics and from quantum fiuctuations. 
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I. INTRODUCTION 

^ ■ Over the last decade, experimental realizations of strongly correlated quantum systems have led to 

■ ' the possibility of studying non-equilibrium quantum processes on a microscopic level. From a theoretical 
point of view, the description of such processes is most challenging because it requires a thorough study 
of the spectrum and correlation functions. 

In this work, we consider a model which describes the hyperfine interaction of an electron spin (the 
central spin) in a quantum dot with a bath of nuclear spins in the dot. The resulting Heisenberg 
exchange interaction is dominant for short time scales up to Ims^ before other mechanisms like spin- 
^ ' orbit coupling or dipole-dipole-interactions between the bath spins set in.^^^ This is an ideal system 
Q \ to generally understand the decoherence of a qubit which is realized by the electron spin,^ and in this 
O . context the loss of quantum information. Many important contributions on this central issue have been 
made by a number of authors using different methods, '''^"^^ as also outlined in the reviews [20,21]. All 
those works rely on often very sophisticated approximation schemes to study the time evolution of the 
^ ' central spin directly. In this work, our goal is to use the exact solution of the model to study its spectrum 
OO , and static correlation functions in the ground state and excited states also in comparison with a simple 
' classical approximation. In the future, this knowledge can be used to obtain exact information about 
non-equilibrium dynamics like the decoherence process. 

The central spin model (or Gaudin modeP^'^^) we consider here describes the isotropic Heisenberg 
On ' coupling of the central electron spin with inhomogeneous exchange couplings to a bath of iV;, 
, nuclear spins Sj=i^,,,^Ni- The nuclei are assumed to be spin-1/2 particles and their couphng gn to the 
external magnetic field h is assumed to be much weaker than that of the electron, ^c, 



H = Y,A,S^- S, - hg^S'^ ~ hg^ ^ 5| . (1) 



X 

b . 

5t 1 The couplings Aj are proportional to the square of the electronic wave function at the positions of the 
nuclei. For a realistic distribution of the Aj, we can think of the index j as measuring the distance 
from the center of the dot. The methods we use in this work, especially the classical approach and the 
integrability, do not depend on the choice of couplings Aj, but for dcfinitcness we assume a harmonic 
trapping potential for the electron. This results in a Gaussian decay of the couplings^° 



Aj = a cxp 



-{jB/Nl'^'f] , (2) 



where the normalization a = xiN^/ Y^^^Li^^V ^iJB/N^^^)'^ is chosen such that the mean value (or 

first moment) xi of the Aj is fixed and the dimension is taken D = 1. Here, the parameter B controls the 
degree of inhomogeneity. We will choose i? = 2, xi = 2 as generic values for inhomogeneous couplings 
and B = 2/5, xi = 2 as an example for nearly homogeneous couplings in numerical diagonalizations in 
later sections. 



We calculate the spectrum, the magnetization profile (5*1) of the nuclear bath spins, and the two-point 
functions {Sq ■ Sj) of the model (1). It is possible to distinguish two types of contributions in these 
quantities: On the one hand, terms appear that can be obtained from a purely classical approach. 
Additionally, we identify terms stemming from quantum fluctuations. Most importantly, classical and 
quantum terms can be of the same order in the two-point function. 

This paper is organized as follows. In Sec. II, we show how to obtain one- and two-point functions 
from the exact Bethe Ansatz solution for the eigenvalues and eigenstates of the Hamiltonian (1). The 
central spin and BCS pairing models arc linked by their integrability which provides a way to calculate 
the magnetization profile (Sj). Two-point functions arc given by derivatives of the energy with respect 
to the Heisenberg coupling constants. 

In the third section, we evaluate one- and two-point correlation functions based on a classical picture 
which, for finite magnetic field, assumes spontaneous symmetry breaking in the model (1), similar to the 
superconducting phase transition in the closely related BCS model. The local magnetization obtained by 
completely diagonalizing the quantum mechanical model with 16 spins agrees very well with the classical 
results. However, for the two-point function, the agreement is less good, which indicates that quantum 
fiuctuations arc of the same order as the classical terms. 

We are thus led to study the exact solution in Sec. IV, especially in order to obtain quantum mechanical 
contributions to correlation functions. This is done for zero and finite magnetic fields. The connection 
with the classical approach is also established. The paper ends with an outlook. 



II. EXACT SOLUTION, LINK TO THE BCS MODEL AND CORRELATION FUNCTIONS 

A. Exact solution 



With a special focus on the magnetic field terms, we rewrite Eq. (1) as 

H = ^^AjSa ■ Sj - HoSq ~ /itS'tot 
hgn and the total polarization S^^^ ~ X)^'' 



(3) 



M. where N ^ Nb + 1 



with IiQ ~ h{gc — (7n), — hg^ and the total polarization S^^^ ~ '^'j=o "^j " T ^ -"^ ' wnere i\ ~ jvt, 
is the total number of spins and M is the number of flipped spins compared to the ferromagnetic all-up 
state. Note that S^^^ commutes with the Hamiltonian, 



(4) 



and thus S^^^. = N/2 — M is a constant of motion. This means that the last term in (3) provides an 
additive constant which we will drop in the following unless otherwise stated. 

The model (3) has been solved by Gaudin^^'^"^ using a coordinate-type Bethe ansatz; an algebraic 
solution has been given by Sklyanin,^^ and is also described in Ref. [25]. The exact solution has been 
used in Ref. [26] to calculate non-equilibrium dynamics in a fully polarized bath. Using the notation from 
Ref. [26] , the eigenvalues A in a sector of given M read 



A 



^ Mb 



NbXi 



k=0 



2 



(5) 



where Xi is the mean value of the Aj and Mf, 
by the Bethe Ansatz equations (BAE) 



M — 1. The set of the ojk, k — 0, . . . , Mi,, is determined 



1 



Nb 

.7 = 1 ^ 



Mb 
k'^k 



^k' 



Uk' - UJk 



2ho 

L^k 



0. 



(6) 



Gaudin^^ showed that there are = N\ / {M\{N — M)\) sets of solutions {wq, ■ • ■ , wa/^} to these equations 
in each sector of given M, one for each eigenvalue A. The corresponding energy eigenstates with a fixed 
number M of flipped spins are given by 



\M) 



Nb 



A, 



nM 



k=Q 



A-'i 



|0), 



(7) 
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where |0) is the fully polarized state | -ff; t, • ■ • , T)> and the arrows ff, -|| for the central spin and t, i for 
the bath spins are used. The normalization factor um was conjectured by Gaudin^'^'^^ and proved by 
Sklyanin^^ for /iq = 0: 



'A/ 



(-1)^^ detX 



Mkk' = k ^ k'. 

[ujk ~ LOk' r 

In Rcf. [26] evidence was given that this holds for finite Hq as well. 

Let us now come back to the eigenvalues in Eq. (5). Due to the HcUmann-Feynman theorem, ^^"^^ two- 
point correlators between the central spin and a bath spin in an eigenstate are obtained as the derivatives 
of the energy eigenvalues 

{So-S,)^dA,A, (8) 

and the expectation value of the central spin polarization is given by 

(S^) - -a,„A. (9) 

By solving the BAE (6) as a function of the couplings Aj and the field ho it is therefore possible to obtain 
the expectation values directly. 

In order to also calculate the magnetization profile (Sj), j — 1, Ni, wc have to use some additional 
features of the integrablc structure of the model (3), as will be described in the remainder of this section. 
Let us rewrite Eq. (3) in the original notation used by Gaudin,^'^ 

Of ■ J-j 



H, = - y] ^i^^h,st, (10) 

j= 

such that wc recover Eq. (3) with ft-t = for 



j=o,j^e 



A, = 1/e, (11) 

and ^ = 0, Eo = in Eq. (10). As pointed out by Gaudin,^^ 

[Hi,He-]^0 (12) 

which means that an integrablc Hamiltonian can be constructed as a linear combination of N mutually 
commuting conserved quantities 

H-.^Y. '^H^ = -'^o E '^^^ - 2 + 2 ^ ^^^^ 

i=o e=o e=o 

with Stot = Z^fco ^e- 

The model (13) has the same eigenstates as the original model (3), even though these are not necessarily 
in the same energetic order. For the local expectation values, one can apply the Hellmann-Feynman 
theorem^^""^^ to the eigenvalues A of H: 

(S^) = -i^a,A=-la,-.A. (14) 
In order to calculate A, we use Gaudin's result^^ for the eigenvalues A*^^^ of in Eq. (10) 

Aw^ly^ i y _J_^!^ (15) 
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with 

Ek = IM, (16) 

such that 



A = ^£,aW 



1=0 

l^Si NtjNt + 1) ho^ 

= 2^.^7:7^ 8 Yl^''- 

e,k " f=o 

By rewriting the BAE (6) in terms of the e^, E-^ and defining g := l//iOj we arrive at 

Ei^-^Ei^^+ro, (18, 

Observing that ejl{Eu - Sj) = Ek/{Ek - Sj) - 1 and Ek/{Ek - Ek') + Ek,/{Ek' - Ek) ^ 1 we can 
ehminate the first term in Eq. (17) by multiplying (f8) with Ek and then performing the sum over k. 
Hence Eq. (17) becomes 

A = ho}_^Ek~ 4 8 2 ' ^ ' 

k=o e=i 

which yields, together with Eq. (14), 

iS]) - 2 ~ ^ = 2 - ^VA. E -• (20) 

fc=0 fc=0 

In summary it is therefore possible to express the two-point function (8) and the local magnetization 
(20) in terms of the BA numbers of the exact solution, which is the main finding of this section. These 
quantities will be analyzed in detail in sections III and IV. 



B. Link to the BCS pairing model 

It is possible to relate spin with fermionic operators, using Anderson spin-1/2 pseudospin operators 

J.... J 



25,33 



S- - - l-c]^c,T-c],c,J, (21) 



S.. = c\A , S+ = CjiCj^ , (22) 



2 

which preserve the SU{2) commutators [5+, 5^] = 2S,jS], [S^,Sf] = ±6ijSf.^^ 

A BCS-like Hamiltonian can be defined by rescaling the integrable model H from Eq. (13) 

Fbcs ■-j^H+-l^e, + — (23) 

where 5*^^^ — Nh/2 — Mh— 1/2 is the conserved quantum number from Eq. (3). In terms of spin operators 
Hbcs therefore reads 



i^Bcs = - E ^7' (^1 -I)- ^(^*°')' + ^"^tot (SLt + 1) • (24) 
Replacing the spin-operators through Eqs. (21) and (22), one arrives at the fermionic representation 

Hbcs = 2 E ^^^L^^*^ ~ ^ E cli^T^jTCji, (25) 

o-=T,J- 
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FIG. 1: (color online) Left panel; Magnetization profile in a quantum dot with Ni, = 15, infinite central magnetic 
field, fto ^ oo, and total field As > 2ht > Ag. The black crosses denote the coupling constants chosen according to 
Eq. (2) with xi = 2, B = 2. Right panel: Corresponding electronic occupation probability according to Eq. (26) 
for free Fermions, with (black crosses) the single-particle levels ej , j = 1, . . . , Nt, according to Eqs. (2) and (11). 
The system is in a two-particle excited state, where the pair occupying the lowest energy level eo = is shifted 
to the energetically lowest state above the Fermi level. 



with the doubly degenerate single-particle levels and the pairing amplitude g ~ Xjh^?^ 

The Hamiltonian (25) describes M pairs of fermionic particles interacting via an attractive pairing 
potential, thus affecting the N doubly degenerate energy levels Cj. In a series of papers, Richardson^^"""^ 
used it to describe pairing in nuclei. In the more recent past, the exact solution of this model has been 
rediscovered to study ultrasmall metallic grains in their superconducting phase. In the thermodynamic 
limit, the solution of the model (23) yields the mean-field BCS solution;^^ we will come back to this point 
in Sec. III. 

From Eqs. (14) and (23) it follows that the occupation probability (rij) := (c^|Cjf -I- c^j^Cjj)/2 of the 
single particle level Ej reads"'"'"'^ 

-1-{SJ), (26) 

which is consistent with Eq. (21). Thus the single particle occupation numbers in the pairing model are 
directly related to the local polarization of nuclear spins in the central spin model. 

We will compute two-point correlation functions and the magnetization profile for different parameter 
regimes in the following sections. For illustrative purposes, let us first check the extreme limits Iiq ^ 0, oo 
in Eq. (26) for the ground state in the sector S^^^. = (this implies that we take N to be even here). In 
the BCS-model, this corresponds to the case of half filling, where the number of electrons 2M equals the 
number of free particle levels TV. For ho 0, the model (3) is S'C/(2)-invariant, so {Sj)\ho=o — 0. Since 
g = 1/ho, the pairing potential is infinitely strong in this limit, such that all levels are occupied and only 
ideal Cooper pairs exist, where each level is occupied by half a pair. 

In the opposite limit, Hq — > oo, the central spin is frozen along the z-direction. The directions of the 
bath spins are simply given by the competition of the antiferromagnetic exchange in Eq. (3) with the 
magnetic field ht- Therefore, all outer bath spins with coupling Aj < 2ht are aligned with the field and 
the central spin, while the inner ones point in the opposite direction. The resulting magnetization profile 
is sketched schematically in the left panel of Fig. 1, where we chose ^(7V6-i-i)/2 > 2/it > A(^j^_|_i)/2+i, such 
that /Sjot = for illustrative purposes. 

For the BCS model (25) this means that the highest states where ej is largest (i.e. Aj = l/sj is 
smallest) are unoccupied. This is the filled Fermi sea for the non-interacting Fermi gas. The level Eq = 



5 



is special in the sense that it is unoccupied in the ground state of the central spin model, which is an 
excited state in terms of the BCS Hamiltonian. From this we conclude that the ground state of the 
central spin model for S^^^. = corresponds to an excited state of the BCS model where the energetically 
lowest pair is shifted to the top of the filled Fermi sea. This is illustrated in the right panel of Fig. 1 and 
will be further discussed in Sec. IV C. 



III. THE SCREENING CLOUD FROM A CLASSICAL POINT OF VIEW 



In this section, we develop a classical picture for the energy and the magnetization profile of the model 
(3) for finite magnetic fields, which turns out to be closely related to the mean-field BCS solution"'^ of 
the pairing Hamiltonian (25). 

It is reasonable to expect that for large coordination number A^f, ^ 1, a classical approach to the 
Hamiltonian (3) yields valuable insights into the physics of the model. The classical approach consists 
of replacing quantum-mechanical spin operators Sj by classical vectors (mj). Especially, for states with 
the same quantum number S^^^, an expectation value (Sj) ^ implies that in this limit, the Hamiltonian 
symmetry (4) is spontaneously broken. This mechanism is analogous to the superconducting phase 
transition in which particle number conservation is broken, (Cj|cj,|) 7^ 0. 



A. Magnetization pattern in the central spin model 

Let us begin by parameterizing each spin in polar coordinates, rrij = 
^{cosifij sin i?j , sin iy9j sin iJj , cos i?^ ) such that |mjp = i for j = 0,...,7Vf,. Our aim is to derive 
the ground state configuration described by the angles tpj, "dj for a given total magnetization S^^^^ and 
fixed central field /iq. 

The classical energy as a function of the azimutal angles Lpj is always minimized by choosing tpo — fj = tt, 
corresponding to antiferromagnetic alignment in the xy-plane. The resulting classical model for the polar 
angles analogous to Eq. (3) is then given by 

^ Nb ^ h 

= i E (^0 + ) - y cos ^0 - y E ' (^7) 

and the total magnetization can be determined from 

2S^,,=J2cos^,. (28) 

The first antiferromagnetic term in Eq. (27) is minimized by large polar angles ?9o + "^j ~ ""j spins 
lying in the xy-plane, while the field tends to keep the polar angles small, analogous to the situation in a 
two-dimensional Heisenberg antiferromagnet with a central impurity.'*^ For finite fields the central spin 
typically acquires a relatively small but finite polar angle, while the bath spins cant into the opposite 
direction out of the plane with polar angles that are closer to 7r/2. Depending on the overall magnetic 
field this results in a characteristic magnetization profile: Those bath spins which are coupled strongly 
are aligned antiferrogmagnetically to the central spin (i.e. against the field), while the more loosely 
bound bath spins at the edge of the dot are aligned ferromagnetically. Depending on the parameters the 
total magnetization Sl^^ is often quite small or even negative. A typical resulting magnetization profile 
is sketched in Fig. 2. 

The minimal values of the angles are most easily found by requiring that the total value of the torque 
I To I experienced by the central spin from the central field and the bath spins has to vanish 

|to| - c)^„i7ci = (29) 
^ ho sin ^0 = 2 E sin(t9j + tJq) - sin i?o ■ (30) 
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FIG. 2: (color online) Schematical orientation of the classical spins according to Eqs. (36) and (37). The fat 
leftmost spin is the center of the dot. The central field leads to a canting of the central spin, which, due to the 
antiferromagnetic exchange, leads to an opposite canting of the neighboring spins. Since an overall magnetic field 
is included which fixes the total magnetization, a non-trivial magnetization profile results. 



Equally well, the torque on each individual bath spin is zero in equilibrium 

|t,| = d.j)^H,i = (31) 
/it sin 7?, = ^sin(t9, +79o) {] ^ 0) ■ (32) 

Obviously, Eqs. (30) and (32) are trivially fulfilled when 'dj=o,...,Nb are multiple integer values of tt. We 
exclude these solutions here, because generally, they do not correspond to minima of the energy, as can 
be seen from the Hesse matrix of second derivatives of Hd . 
We now insert Eqs. (30) and (32) into Eq. (27) and obtain 

From Eqs. (30) and (32) it follows that if the fields /io,t are given, then we can solve for the angles "dj, 
which are given by 

tan^, = -^, (34) 

where 

(5 = tan?9o5 v — 2ht/ cos-Oq. (35) 

The angles in Eq. (34) shows the generic behavior described above unless extreme values of the parameters 
arc assumed: The magnetization changes from alignment with the field for the outermost bath spins 
{Aj — > 0) through the xy-plane {Aj ^ v) to near antiferromagnetic alignment for the most strongly 
coupled spins near the center [Aj > v). 

The components of the magnetization along the field and in the plane can be found explicitly 
by using 

tani^o = m^/m^ = S (36) 
S A 

tani^, = mf/m^. = ^— for j = 1, . . . , Nb, (37) 

v ~ Aj 



from which it follows that 



1 J. 

mn 



2^{v~A,Y + {A,8f' 2^(v^A.i)-^^{A,5f 



(38) 
(39) 



Similar equations were obtained using methods of classical intcgrability in Ref. [44] . In order to determine 
the parameters 8 and v we obtain from Eqs. (38) and (39) for the total magnetization along the field 



1 A 

1 V — Ai 



2S' = , + V ' =2N - M . (40) 
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Equation (30) for the central field now reads 



Eqs. (40) and (41) fix 5 and v uniquely for a given S^^^ and hg, so that all classical vectors arc known, 
which is the central result of this section. 

Finally, one obtains the corresponding expression for the energy from Eq. (33) without the trivial 
/it-term 



1 



iu-A,y + iA,sy 



1/2 



A' . (42) 



This parametrization of the ground state energy in terms of v, 6 and the Aj will be helpful in separating 
classical from pure quantum contributions in the exact solution later on in Sec. IV C. 

It is interesting to note that an alternative derivation of Eq. (41) is obtained by considering the 
magnetic fields /ig.t as canonically conjugate to mg, 1 so that ho = 9mg X]j=i '^o ■ nT-j and ht = 

The classical spin-spin correlation function between the electron and nuclear spins can be obtained 
from Eqs. (38) and (39), namely 

(l + 6^)A,-i^ , , 

mo-m, = ^ , : =■ (43) 

B. Connection with tlie BCS-model 

Very similar relations were derived^^''**''^^ for the thermodynamic limit of the BCS-pairing model (25) 

A2 N, 

9 1 

- = y , (45) 



N-2M = y , ^" ^ - . (46) 

Here A is the superconducting gap, /i the chemical potential and iJgcs ^^^^ ground state energy of (25) 
in the thermodynamic limit. Eqs. (45) and (46) are equivalent to (40), (41), if Eq. (11) and the following 
relations hold, 

6^A/^l, v=l/pi, ho = l/g, (47) 

and if furthermore, the sign of the j = 0-term in Eqs. (45) and (46) is changed. The latter condition 
reflects the fact that the ground state of the central spin model corresponds to a special single-pair excited 
state of the BCS-model. This point will be discussed quantitatively in Sec. IV. 

The mechanism of spontaneous symmetry breaking in the classical/mean field approach is completely 
equivalent in both the BCS and central spin models. In order to see this, we use the pseudospin repre- 
sentation (21), (22) to write the BCS gap in the pair-excited state corresponding to the ground state of 
the central spin model as 

A = gymf. (48) 

Inserting the last of relations (39) and substituting Eq. (11) and the first two equations from (47), one 
re-obtains the gap equation (45). 
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An important difference to the BCS-solution consists in the order of magnitude of ho = 1/g. To 
obtain a well-defined energy per particle in the thermodynamic limit, 1/g = 0{N) scales with the 
number of particles. In the quantum dot. however, the experimental situation corresponds to Hq = 0(1), 
thus not scaling with any extensive parameter. It is instructive though to consider the limit of infinite 
central magnetic field, shown in Fig. 1. Then 6—^0 and from Eq. (41), ho w J2f=i |2(£j — ^ 1- 

Furthermore, Eq. (40) yields Ajv/2 < < Afff2+i, so that Eqs. (38) and (39) reproduce the magnetization 
profile shown in Fig. 1. In this extreme limit, quantum fluctuations are suppressed completely and the 
classical picture is exact. Accordingly, the classical mean field approximation is generally better justified 
for the BCS model. However, for general fields kg = 0(1), apart from the classical contribution discussed 
in this section, important quantum fiuctuations will occur as well as will be shown in the next section. 



C. Analytical results: Small field limit 

Eqs. (40) and (41) can be solved numerically to determine the parameters i^, S from which the mag- 
netization profile Eq. (39) and the two-point-function Eq. (43) are obtained. However, in the physically 
most relevant limit of small central fields and large particle numbers it is useful to derive approximate 
analytical expressions for the one- and two-point correlators. Therefore, we will first calculate the pa- 
rameters S, V from Eqs. (40) and (41) to leading order in ho, before inserting these results into Eqs. (39) 
and (43) for the correlation functions. 

According to Eq. (41), a small central magnetic field corresponds to 

where the index 5^^^ is the leading term of (5 in a small-field expansion of 5. Since we derived Eqs. (40) 
and (41) for a large number of nuclei, we restrict ourselves to the terms leading in Nb here. Eq. (49) is 
consistent with Eq. (47): Both imply that a diverging pairing strength in the BCS-pairing model leads 
to a diverging superconducting gap. 
In the same limit, Eq. (40) leads to 

25<'^^Sl,^N -vNbX_^, (50) 

where we defined the moments X£ 

Nbxr.= Y.^'r (51) 

The moments with negative (positive) integers i are determined predominately by the smallest (largest) 
coupling constants. 

We consider here a sample which is not macroscopically polarized, i.e. S^^^ = 0(1). The case of 
macroscopic polarization will be dealt with in Sec. IV. Together with Eq. (49), we then obtain for the 
leading term of u for small fields, v^^^: 

Making the same approximations in the expression for the classical ground state energy, Eq. (42), and 
inserting Eqs. (49) and (52), we obtain the leading term for small ho 

For small central fields, this yields the following expressions for the leading terms in a large- TV-expansion 
of classical one- and two-point correlation functions in the ground state: 
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FIG. 3: (color online) Magnetization profiles as a function of the site j and the central field ho for S^^t = with 
the couplings (2) where xi = 2, B = 2, obtained by inserting the numerical solution of Eqs. (41) and (40) for 5, v 
into Eqs. (38) and (39) for m|. The short blue lines denote the leading contribution for small fields and large 
particle number, Eqs. (55) and (56). 

where is the leading term in an asymptotic expansion of in the inverse particle number, := 
\IA(xN\^Ax^ and A(xiVb) = Aj is treated as a continuous function of x. 

D. Quantitative comparison with numerical results 

We illustrate the classical results in Fig. 3, where magnetization profiles n^j^o are shown, after 
solving Eqs. (40) and (41) numerically for S^^^. = and different fto- The small-field asymptotes from 
Eqs. (55) and (56) are depicted as well. 

We now discuss the question to what extent these classical expressions can be identified with the 
quantum-mechanical expectation values for large particle number and small central field. 

In order to do so, we first compare our results with a complete diagonalization study for a system 
with TV = 16 sites as an additional independent check. The coupling constants in this system were 
chosen according to Eq. (2) with xi = 2, B = 2. In the next section, we will see that the complete 
diagonalization study also enables us to classify low-lying excited states according to the distribution of 
the corresponding BA roots, which is not possible a priori. 

In Fig. 4 we compare the diagonalization results with the full classical expressions Eqs. (38) and (39) 
and with the approximate results (55), (56) for three different values of Hq. The small- field expansion 
(56) deviates from the exact data essentially at large distances from the center of the dot, where the more 
weakly bound spins are located. On the other hand, Eq. (39) with values for 6, v obtained by solving 
Eqs. (40) and (41) numerically deviates from the exact solution only by a few percent or less. 

Comparing the classical expression for the two-point function Eq. (43) with the diagonalization results, 
one notices considerable differences, see Fig. 5. 



10 




5 10 

3 



FIG. 4: (color online) The local magnetization {Sj)o obtained for a. N = 16-spin system with couplings as in 
Fig. 1. The central field assumes values ho = 0.5, 1, and 2, as indicated in the figure. The total polarization is 
fixed at Sfot ~ 0. Data from complete diagonalization (red crosses) are compared to the small-field expressions 
Eqs. (55), (56) (green diamonds) and the mean-field result Eq. (39) (blue circles), where S, v were obtained by 
numerically solving Eqs. (40) and (41). 
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FIG. 5: (color online) Two point function {So ■ Sj)o in the ground state for the model with the same exchange 
couplings chosen as in Fig. 1 and central field values ho = 0, 0.5, 1.0, 2.0 from below. Black crosses were computed 
from a complete diagonalization. Green squares were obtained from Eq. (93). Red diamonds stem from the 
classical solution, Eq. (43). Blue circles are based on Eq. (110), where classical and quantum contributions are 
summed. 
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IV. CORRELATION FUNCTIONS FROM THE EXACT QUANTUM-MECHANICAL 

SOLUTION 



Whereas in the previous section a classical picture of the central spin model was sketched, this section 
contains a systematic study of the exact quantum-mechanical solution, where the contribution of quantum 
fluctuations to the correlation functions will be emphasized. We will first obtain approximate analytical 
expressions for correlation functions in the regimes of zero and weak central magnetic fields, before 
recovering the classical picture from the previous section in the appropriate limit. 



A. No field 



For vanishing magnetic field, the Hamiltonian (1) is S'[/(2)-invariant and commutes with all components 
of the total spin, [H\h=o, S'tot] = 0. In other words, all states within one spin multiplct, obtained by acting 
with S'jQt on highest weight states, are energetically degenerate. In the expression for the eigenstates (7), 
application of S^^^ corresponds to choosing ujk,u = 0. Indeed, for ho = 0, it is easy to see that if 
{wi, . . . , Wm} is a solution of the coupled set of equations (6), then {wi, . . . , lja/, 0} is a solution as well. 
Both these solutions are energetically degenerate, according to Eq. (5). This situation is analogous to 
the Heisenberg chain, where sets of only finite roots encode the highest weight states. "'^■^^ Here and in 
the following, the ground state energy and expectation values in the ground state will be labeled by the 
subscript o- 



1. Ground state 

The ground state maximizes 'Y^^lo^k- It turns out that the corresponding highest weight state has 
Mb = 0, so that only one Bethe number has to be determined from 1 — Y^^li uiq-a- ~ ^- ^"^^ 
interested in the energy levels for large particle numbers. In the ground state, wq = 0{N), which allows 
to rewrite the single BA equation in terms of the moments Xn defined in Eq. (51), 



oo 



A^^ES = 0- (57) 



' 1 '^o 

n— 1 ^ 



We define i^oNb, such that both Xn and llIq are 0(1). Then, according to (5), the ground state 

energy reads 

i^o = -^f^o-:^) (58) 



2 V " 2 

and successive orders of ujq in an asymptotic expansion for large particle numbers can be obtained by 
inverting Eq. (57) order by order. Let cDq"'' be the expansion of ujq in powers of N^^^ up to order n, i.e. 

limAT^^oo (^Q — '^o"'') N^'^^ = 0(1). For n = 3, with d := {NhXi)~^, ?/„ :— Xn+i/xi, wc obtain 

~(3) 

^ = l + yid+{-yl + y2)d^ + {2yf-3yiy2 + y3)d\ (59) 
x± 

This leads us to conjecture that the coefficient of d" in the expansion of {ojq/xi — 1) is given by the 
nth coefficient in a Taylor expansion of nl In [(l)(d)/xi] in the variable d, where the generating function is 



Before continuing, let us make two comments on Eq. (59): To begin with the leading term, when plugged 
into Eq. (58), yields the overall ground state energy of the classical model (27) with no fields, where the 
central spin is pointing in the direction opposite to the nuclear bath spins. Finite-size corrections, given 
by the sub-leading terms in Eq. (59) therefore represent quantum effects. Secondly, for the homogeneous 
model Aj = AVj, all but the first two terms on the right hand side of Eq. (59) vanish. 

It is now straight-forward to evaluate the moments for a given distribution of the coupling constants 
Aj. For the particular case of the choice in Eq. (2) it is possible to use the Euler MacLaurin summation 
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FIG. 6: (color online) The two-point function (So ■ Sj) for zero magnetic field in the ground state (gs) and the 
lowest ten excited levels (couplings as in Fig. 4), obtained from complete diagonalization. The circles in the left 
panel are the analytical result Eq. (64) for the ground state. 



formula to find an expansion of the moments in the parameter d. Writing yi ~ y'f'^ + dy^i \ we find 



(0) 

Vi 



2xi 



B Erf ( V2 B) 

2xlB^ , 4^2 Erf(^/2S) 



7rErf^(2) 



7rV2Erf^(B) 



(60) 
(61) 



Therefore, the coefficient of in Eq. (59) becomes 



(0) 



(0) 

V2 



2x1 
7rErf^(2) 



-Erf2(\/2S)/Erf(B) + 2Erf (Vs S)/\/3 



(62) 



As expected, this latter expression tends to zero for B ^ 0, which is the homogeneous limit in the 
couplings (2). 

From Eq. (58), one then obtains for the ground state energy 



(0) 



,,(0) 
2/2 



(63) 



The first term is the classical result, where the central spin is aligned antiferromagnetically with respect 
to the bath spins. The second and third terms constitute quantum corrections. 

In order to calculate the two-point correlation function {Sq ■ Sj)o in the ground state, one combines 
Eq. (8) with Eq. (63). This yields up to order d 



dA, 



(64) 



Again, the leading contribution refiects the classical picture of antiferromagnetically aligned spins. For 
zero central field, quantum fluctuations lead to a non-trivial dependence on the distance j between the 
bath and the central spins. This is a pure quantum effect, as can be seen by comparison with Eq. (54). 

Especially, quantum fluctuations decrease {Sq ■ Sj)o below the classical result —1/4 if Aj > yf^ /2, i.e., 
for the strongly coupled bath spins. In Fig. 6 the result (64) is compared to complete diagonalization 
data. In the homogeneous case Aj ~ AVj, Eq. (64) reduces to the result found in Refs. [48,49]. 
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FIG. 7; (color online) The local magnetization {Sj) for A^;, = 15 bath spins, with the couplings (2) where Xi = 2, 
B = 2. Black crosses stem from a complete diagonalization, blue circles from Eq. (67). For the central spin, 
(5*0) was determined such that Sfot is fixed, which in leading order is given by Eq. (68). From top to bottom, 
= 7,6,5,4,1. 



To calculate the magnetization profile (5'J)o, we can employ Eq. (20) together with Eq. (59). Let us 
first consider the case of a fully polarized bath. For S^^^ = N/2 — 1, in leading order this leads to 

(5|)o = l-d'Al (65) 

{S^)o = -l + dy^^\ (66) 

At smaller values for 5^^^, we can still use Eq. (20), keeping in mind that it has been derived at finite 
Hq. Wc thus have to perform the derivative in Eq. (20) before taking the limit ho 0. In this limit, the 
eigenvalues A in Eq. (19) were given in Ref. [50], Eq. (39). Using that result we obtain the expectation 
values in the respective ground state of each sector S^^^ 

{S^)o = ^ (1 - Af) , j^l,...,N, (67) 

(^o)o = ^(-l + 2dy(")). (68) 

For ^tot = A^/2 — 1, Eqs. (65) and (66) are recovered. In the opposite limit, Sq ~ 0, the polarization 
vanishes, as expected from the S'[/(2)-invariance in this case. In Fig. 7, we compare the analytical results 
(67), (68) with complete diagonalization data. These illustrate the fact that the magnetization profiles 
are different for energetically degenerate states. 



2. Excited states 



In the classical picture, the lowest excitations at Hq ~ above the Nccl-likc ground state are created 
by flipping the spins in the outer region where the exchange with the central spin is weakest. In the exact 
solution of the quantum mechanical problem (5), (6), excitations can be of two types: 

• Spin excitations with a change of M , i.e. the number of roots. 
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• Particle-hole excitations, where the location of roots is changed with respect to the ground state, 
but the number of roots is kept fixed. 

For /iQ = 0, both types of excitations are energetically equivalent: Adding a root 7^ is equivalent to 
moving a root ujk = to a finite value. 

Let us consider the excited state where Mb + 1 roots are different from zero. We focus on low-lying 
excitations here, so Mf, does not scale with the particle number. Thus compared to the ground state for 
25*^0^. = N — 2{Mb + 1), there are now Mb additional roots away from the origin. In the set of equations 
(6), there is one root which scales like the particle number; we denote it by luq, i.e. ujq — 0{N). For the 
other Mb roots, ujk = 0{1). 

We define the moments of the additional non-zero roots as 7„ :— X]fc=i ^1! • Performing an expansion 
analogous to Eq. (57), one obtains for the root ujq = 0{N) the equation 



" ,,,n 



(69) 



which again can be inverted order by order. Including terms of order 0{l/Nb), 

ojQ = - - 2-fi + yi + 2dj/i7i - 2^72 + d{y2 - y\)- 
This leads to an expression for the energy in terms of the 7n 



Vi 



(0) 



272 



(0) 

yi 



(70) 



(71) 



Let us look at the simplest case, Mt, 



1. The corresponding equation for the additional root loi reads 
A, 







(72) 







By sketching the Ihs of this equation, one sees that ui is located between two couplings. Indeed, 
for the lowest excitation, we can set cji = A^^ + S^,^. In leading order, we then obtain Sn^ = 
ANb/{J2f=i^ Aj/(Aj — AjvJ — 1) = 0{1/Nb) and J^r^ > 0. One can generalize this result to uji = A^ + Sg, 
as long as Si ~ 0{\/Nb) and 5i > 0, i.e. for £^ I. Then Eq. (64) is modified according to 



{So ■ S 



_ J -\ + ly^^^-dA„ j^. 



(73) 



This corresponds to the classical picture of spin flips with respect to the ground state at the outer edges 
of the quantum dot. That result generalizes further to the case of more than one excitation, Mb > 1. If 
more than one BA root is present, different root patterns are possible. Let us call the distance Aj — Aj+i 
the jth coupling interval. We call an interval occupied if one root is located within this interval. 

One type of root configurations consists in only real roots and occupied intervals, with no consecutive 
occupied intervals. Another type of root configurations involves consecutive occupied intervals. However, 
depending on the special choice of the coupling constants and the central magnetic field, roots in such a 
configuration can be driven into the complex plane, thus forming complex conjugate pairs. '^^'^^'^^ 

From these observations we conclude that the two-point function {So-Sj)cx yields significant insight into 
the underlying root configuration of a low-lying excited state. Vice versa, if the root configuration for low- 
lying excitations is known, the corresponding two-point function can be predicted at least qualitatively. 
This prediction confirms nicely the physical expectation. 

In Fig. 6, we depict {Sq ■ Sj)cx for the lowest nine excited states for — 16 particles with the 
couplings chosen according to Eq. (2) with xi = 2, i? = 2. The data have been obtained from complete 
diagonalization. The analytical result (64) for the ground state is given as well, from which the analytical 
predictions for excited states are obtained straightforwardly by changing the sign of the corresponding 
spins, like in Eq. (73). 

It is instructive to consider the corresponding root configurations of those lowest nine excited levels. 
These are shown in Fig. 8 for the highest weight states, i.e. without roots in the origin. The physical in- 
terpretation of the root locations as spin flips with respect to the ground state is revealed when comparing 
the root pattern level by level with the j-dcpendcncc of the two-point function. 
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FIG. 8: (color online) The location of the additional BA numbers which are 0(1) (red crosses) for the lowest nine 
excitations (top left: level 1, bottom left; level 5) with Nb = 1^, xi = 2, B = 2 in Eq. (2) for the couplings. The 
six smallest couplings are shown here (blue dots). Not shown is wo = 0{Nb). 



This interpretation carries over to the magnetization profile. We shovif those magnetization profiles 
corresponding to the lowest nine excited levels in Fig. 9. However, as in the ground state, an important 
diS'erence consists in the degeneracy of (^o • Sj)cx for all states within one multiplet. Whereas the 
two-point function is independent of the total magnetization S^^^., the local magnetization (5|)ex does 
depend on that quantity. In Fig. 9, we only give the magnetization profiles for the highest weight states 
parametrized by the roots sketched in Fig. 8. By adding additional roots in the origin, i.e. by lowering 
iSfot, the two-point function is not altered, but (S'|)cx is changed by an overall prefactor like in Eqs. (67) 
and (68). Namely, proceeding similarly as in the derivation of Eqs. (67) and (68), one obtains the leading 
terms of the magnetization profile in low-lying excited states 



(5| 



N-2n 

N-2r, 



1 



1 



- 2(PA^) 



(74) 



Here n > 1 is the number of non-zero roots. For low-lying states, these are located close to the couplings 
Ae-^ , ■ • . , Ag^ , as can be seen from Fig. 8. All states with the same n but different S^^^. are energetically 
degenerate, i.e. have identical two-point functions {So ■ Sj)cx, but different magnetization profiles (S'|)ex- 
For the highest weight states we have S^^^ = N/2 — n, which is also the case for the root configurations 
depicted in Fig. 8. According to Eq. (74) the magnetization profile therefore can be read off from the 
number and location of Bethe roots. This is confirmed by complete diagonalization data shown in Fig. 9. 
Deviations from Eq. (74) are due to interactions between the excitations, which were neglected in the 
derivation of Eq. (74). 



B. Weak field 



A finite magnetic field couples to both the electronic and nuclear spins. As described below Eq. (3), 
the coupling to the nuclear spins is trivial and can be accounted for at the end of the calculation. Let 
us thus first focus on ho > 0, ht = 0. As stated in the previous section, for ho = 0, the eigenvalues of a 
given multiplet are degenerate, which in the root pattern is encoded by roots in the origin. For ho 0, 
this degeneracy is lifted due to the broken SU{2) invariance. Thus it is reasonable to assume that the 
zero roots are driven away from the origin by a finite magnetic field ho- In the weak field limit, this is 
confirmed by the large-(7-expansion of the Bethe roots^° of the BCS-Hamiltonian (25), which are related 
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FIG. 9: (color online) The magnetization profile (Sj )gx in the ground state (gs) and the lowest nine excited levels 
corresponding to the root configurations of Fig. 8, obtained from complete diagonalization. The circles for the 
ground state in the left panel are the analytical results from Eq. (67). The analogous results for excited states 
agree with Eq. (74) to the same order. 

to the Bethe roots of the central spin model via Eq. (16). In this section, we will derive the energy 
eigenvalues and expressions for the screening cloud and magnetization profile for small but finite central 
field. 



1. Ground state 

Starting from Eq. (6), we include a finite field ho in Eq. (69). We then arrive at the following equations 
for ujk^o 



OO OO (-> , 



71=1 



A 



— Aj - LUk 



n=0 



Ah 



2E — 



2/i 







= 0. 



= 0. 



k' = 1 



(75) 
(76) 



Thus the only effect of in Eq. (75) compared to Eq. (69) is to add a term — 2/io to Ni,xi, i.e. up to 
order 0{1/Ni,) we have, 

uo = NbXi - 2/io - 271(1 - rfj/i) + yi(l + 2hod) - 272^ -f d{y2 -yl). (77) 

One then obtains for the ground state energy an expression which still involves the Mjy non-zero roots 

(78) 



1 /iQ 71 J/l d r 9 , 1 

Eq = + y + y- y- 2L^2~^i+ "^^ + ^^^^ ^ 



Let us now multiply Eq. (76) by cuk and sum all terms fc = 1, . . . , A/f,. We assume that ho is sufficiently 
small so that max{|wfc|} < ^jv^ and find 



J=2 ^0 



(79) 
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with xo = 1. Here, we aim at calculating the energy up to ©(/iq). In analogy to Ref. 
make the Ansatz 

71 = ciho + C2hl + 0{hl) 

72 = d,hl + Oihl). 



we therefore 

(80) 
(81) 



Then, including terms ©(/iq), the coefficients ci.2 are found by inserting that Ansatz into Eq. (79), 



Cl 



C2 



2Mb 

Nb-1 



x-i 



di+0{(f), 



(82) 
(83) 



where we only keep the leading finite-size terms. 

An additional equation is thus needed to determine di. This is obtained by adapting the techniques 
used in Ref. [50] to our problem. We then find that in leading order in /ig, the roots uJk=i^...^Mb ^-re related 
to the zeros of associated Legendre polynomials. 



L 



Ah 



2ho 



0. 



(84) 



This is a polynomial of degree Mb, i.e. Lj^^''{x) = cYYl'l^{x—2hQ/ujk), where the constant c is determined 
by the asymptotes. Consequently, the logarithmic derivative at a; = is 



IniMfH^;) 



x=0 



71 



(85) 



On the other hand, 



= Mb/{Nb — 1), which in combination with Eq. (85) con- 



firms Eq. (82). Analogously, the second logarithmic derivative \\nL j^^^'' [x) 
Nb)l {{Nb - 2){Nb ~ 1)2), which leads to 

MbjMb + 1-Nb) ^ 72 
{Nb-2)iNb-iy (2/io)2- 

Combining this equation with Eq. (81), one finds in terms of S^^^ 

iNb-iy-Hs,^,y 



Mb{Mb + 1 



di - 



C2 



(iVfc-2)(iVfc-l)2 



-X-i- 



{Nb~ir~A{s,i,y 



Nb, 



{Nb - 2){Nb - 1)^ 

where the latter relation follows from Eq. (83). Then the ground state energy reads 



(86) 

(87) 
(88) 



1 




qz Qz 
^tot o tot 


Ad ^ 




Nb- 




1 


(0) 

Vi 


d 


(1) f (o)\2 


Id ^ 


2 


^ 2 


Vi - [Vi ) 



jdyiho + - -{y2 - vl) 



(0) 
2/2 



Stot/io(l - 2dyi) 



+ 0{d^) 

l-4«ot)^(0).2 



2Nb 



o{d^),{m) 



where we have defined the total magnetization density s^^j :~ S^^^/{Nb — 1). In the last equation, the 
leading orders in a finite-size and small-/io expansion are given. For 5*^^^^ = 0, the central magnetic field 
does not enter linearly, but due to second-order spin-exchange processes only quadratically. 

It is interesting to note that when one adds the additional total field term —htS^^^ from Eq. (3) to 
Eq. (89), the lowest levels E{){Sl^^, ho) display a pattern which is strongly reminiscent of light rays forming 
a caustic in optics. This is visualized in Fig. 10. Thus for a finite total magnetic field ht, a small but 
finite range of values for ho exists where the state with S^^^. = is the ground state. Or, coming back to 
the original Hamiltonian (1), this means that for any finite ratio g-a/go, one can adjust the field such that 
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FIG. 10: (color online) Lowest levels in the sectors M — 7 ... — 7 for Nt — 15 and couplings as in Eq. (2) with 
xi = B = 2, obtained from a numerical solution of Eq. (6). Black lines: ht — (the fat black line is the energy 
for Stot ~ 0); blue lines: ht = 0.035. The inset is a zoom into the region around the caustic. The red arrows mark 
ho values such that the state with ^tot = is the ground state. 



the ground state has a given magnetization S^^^. For the example shown in Fig. 10, the ground state has 
zero magnetization for ht = 0.035, ho « —0.61, which corresponds to gn/Sc ~ 0.054 with hg^ = ht. 

Analytically, the relation between ho and ht for a given Stot is found from ht ~ ^ds^^^Eo, with Eq 
given in Eq. (89). This leads to 

^ht = (1 - 2dy,) + !£m^ ( Jl^y , (90) 



up to higher order corrections. For S^^^ = 0, this is inverted to 

-ho = {Nb - 1)(1 + 2dyi)ht = UVb + - 1 Ut + 0{d), (91) 



which for the numerical values chosen in Fig. 10 yields ho ~ —0.60, in good agreement with the exact 
numerical data from the Bcthc Ansatz. From the classical Hamiltonian (53), only the leading contribution 
in the particle number in Eq. (91) is recovered. From Eq. (91), the ratio of g-factors in the St^t = 0-sector 
is deduced (setting hgn = /it). 



g-n X 



-=Nb + ^ + 0{d). (92) 



This means that when the ratio of electronic to nuclear ^-factors equals the number of nuclear bath spins, 
then an overall magnetic field drives the system into the non-degenerate S^^^^ = state. Since this ratio 
is of the order 0(10'^), it is very realistic to probe this regime in an experimental setup. 

Let us now consider correlation functions. From Eq. (8), we obtain for large particle numbers 

1 d d^ r*°^ /?2 

{So ■ S,)o = -- + -VI dA, - (-2^2 + ivl + iA] - 4yiA,) + ^ ^ + 0{d^). (93) 

with C2°'' = (1 — 4 (stot)^)/^b • If 01^6 is interested in the coefficients of the asymptotic l/A^6-expansion 
one should again apply the Euler-MacLaurin formula as in Eq. (63). We do not want to dwell into these 
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technical but straightforward details here but rather compare the analytical prediction with exact results 
from complete diagonalization. Such a comparison is shown in Fig. 5 in Sec. HID. 

As expected, the approximation (93) is reliable quantitatively only for small fields, according to our 
weak-field assumption max{|a;fe|} < Aatj^. From Eq. (80), one estimates \LUk \ ~ 2ft,o/A^6, which means that 
Eqs. (89) and (93) are valid for \ho\ < An^Ni,/2. For our choice of parameters Ni, = 15, xi = B = 2 in 
the complete diagonalization of Eq. (2), this means |/io| ^ 0.7. But even for larger values of Hq, Eq. (93) 
is qualitatively correct: A finite central magnetic field leads to an enhanced ferromagnetic correlation 
between the central spin and the rather loosely bound bath spins at larger distances from the center of 
the dot, and to an enhanced antiferromagnetic correlation between the central spin and the bath spins 
closer to the center of the dot, which is also consistent with the classical magnetization profile in Fig. 3. 

We also want to compare Eq. (93) with Eq. (54), obtained within the classical picture for S^^^ = 
for ho ^ 0. In that approximation, the equations differ from each other by field-independent terms 
proportional to d. This is understandable: We have seen in Eq. (64) that these terms constitute finite- 
size corrections which stem from quantum fluctuations and are thus not present within the classical 
approach. These lead to an increase of the amplitude of two-point functions. Especially, for the stronger 
couplings, values smaller than the classical bound —1/4 arc reached, a clear sign of entanglement and 
non-commutativity of the quantum spin operators. 

In order to determine quantum fluctuations to (Sj) in the small-field limit, one has to solve the set 
of equations (18) for the Ej in order to determine the eigenvalue in Eq. (19). From this, the local 
magnetization is obtained via Eq. (20). The small-field expansion of the set of Eqs. (18) has been studied 
in detail in Ref. [50]. From that work, it follows that for the ground state of the central spin model at 
fixed Nb, Mb, 

J2Ek= E^, iMb,Nb-l, x^p - 2EP/Nb) (94) 
fe=i 



with 

Eg,{Mb,Nb-l,x^p) = 
From Eq. (77) one computes 



Mb{Nb-Mb) , MbNbX^ 



2ho Nb - 1 

Vbii4(jVb - 

{Nb-lf{Nb-2) 



Oh ( 2 , NbMbiNb~l~Mb) , .q.^ 



Eo = d(l + 2dho + 2dji + dyi) + 0{hi,d^), (96) 

and 7i follows from combining Eqs. (80) and (82). In linear order in /iq and including orders 0{d^), one 
then obtains the magnetization profile 

(5|) = sl,,{l-2A]d^-l6hoA]dHl,, + A]d^{iy,~2A,)) 

+ [A-^ - x_i(l - 2d^A]) + 2d\-i(l - 2d^A^^)) ^~j^^'}°f ho (97) 



with SjQj :~ S^Q^/{Nb — 1). The polarization of the central spin, (Sq), is fixed by the sum rule S, 



^^q(5|). Comparing Eq. (97) with Eqs. (55) and (56), one again recognizes the effect of quantum 
fluctuations which are now sub-leading with respect to the classical contributions. As expected, these 
reduce the amplitude of the magnetization profile, signaling the effects of entanglement. 



tot 



2. Excited states 



In Sec. IV A 2, we found excitations for ho = 0. We can proceed similarly for ho ^ 0. The expression 
(78) for the energy is still valid, but the 71^2 are different now. Let us first consider single-particle 
excitations, parametrized by a single root toi located on the real axis between Aj^i and Aj. Instead of 
Eq. (76) for /c = 1, . . . , Mb, the corresponding set of equations now reads 

l + ET^-2Ef-) -2E^ + — = (98) 
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FIG. 11: (color online) The two-point function {So ■ Sj)cx for ha = 1.0 for the ground state and the lowest nine 
excited states, obtained from complete diagonalization with Ni, — 15, and the couplings according to Eq. (2) with 
xi = B = 2. 



OO 

j=0 j=0 ^ ' ^ j=0 ^ ' k'=2 

A: = 2, . . . , Mfc. 

By multiplying the latter equation with ujk and taking the sum fc = 2, . . . , Mb, one arrives at an equation 
similar to (79), with Mb ^ Mb - I, Nb Nb ~ 2 and a;_„ .t_„ - 2/{NbLj]^), n = 0, 1, . . .. Thus the 
coefficients in Eqs. (80) and (81) are now 

cT - (100) 
= f 2d + A _ TV.^c-i) (101) 

with = dx(Nb-^ ^h~2, Sl^^ 5'fot) as defined in Eq. (87). 

A small field does not change the root pattern of the lowest excited states qualitatively. Very similar 
to the discussion after Eq. (72), one can still make the Ansatz uje = A(^ + 5e for single particle excitations. 
Then for small ho and £^ 1, one again finds that 5e = 0{1/Nb). This picture carries over to multiparticle 
excitations, except that a finite field Hq can lead to complex conjugate pairs of roots. "^^''^"'^'^^ Thus the 
low-energy excitations are still given by approximately independent spin flips of the outer bath spins. 
This is best seen when comparing the correlation functions {Sq ■ Sj)cx, {Sq)ck with the root patterns 
corresponding to the excited states. In Fig. 11, the two-point function {Sq ■ Sj)cx is shown for the ground 
state and the lowest nine levels with a central field ho = 1 in the sector S^^^. = 0. The corresponding 
magnetization profile, (S'J)ox is sketched in Fig. 12, and Fig. 13 shows the underlying root patterns. 

Although qualitatively, the results are similar to those shown in Figs. 4 and 8, there arc two important 
differences. Firstly, the degeneracy between states within one multiplet is lifted, so that both the onc- 
and two-point functions depend on the total magnetization (in Figs. 11 and 12, we have chosen S^^^ = 0). 
Secondly, the ordering of root configurations according to their energies is different. For example, the 
third excited level for /iq = is given by a two-particle excitation (two flipped spins), as shown in Fig. 8, 
whereas for ho = 1.0, such a configuration yields the fourth excited level, cf. Fig. 13. 



„v^/^fc\ „ ^k' 2ho 
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FIG. 12: (color online) The magnetization profile for ho = 1.0 for the ground state and the lowest nine excited 
states, obtained from complete diagonalization with A'^i, = 15, and the couplings according to Eq. (2) with 
xi = B = 2. 



C. From the exact solution to the classical picture 



In this section, we want to make contact with the classical picture presented in section III, starting 
from the exact solution for large N and small polarization, i.e. M = 0{N), and finite central field. In 
this situation, the question arises whether the Bethe roots form a dense distribution in the complex 
plane, which would permit a continuum description. In Fig. 14, we show both the roots tok and the 
inverse numbers iS^ = l/ujk for iV = 16, M = 8, i.e. S^^^ = 0, parametrized by ho- One can show^° 
that for = 0{N), the distribution of the Ek can be described by a cut in the complex plane in the 
thermodynamic limit M, N ^ oo, M/N fixed. However, for the central spin model, we are interested in 
= ho = 0(1). Fig. 14 suggests that such a continuum description still is possible in this case. To see 
this, we first review Richardson's^^ line of arguments for g^^ = 0{N). 

Consider the function 

Mb Nt 

where Ej, Ek are related to Aj, cuk according to Eqs. (11) and (16) and Eq = 0^. Then for g^^ = 0{N), 
the function F is expanded as F = Fq -I- Fi -I- . . ., where F^ = O (A^~''+^) . Such an expansion is justified 
rigourously by showing that F obeys a differential equation which can be solved order by order. We first 
consider the case where in the leading order, all Ek merge to form a cut in the complex plane along an arc 
which is symmetric with respect to the real axis. The endpoints of the arc, a and a*, are parametrized 
by two real quantities, a = fi -\- iA, where fi and A are the chemical potential and the superconducting 
gap of the BCS model (25). This situation corresponds to the ground state of the BCS model. Then 

The ground state of the central spin model corresponds to a particle-hole excited state of the BCS-model, 
where one root Eo is taken away from the arc and instead is located on the positive real axis, close to 
£o = 0~, namely Eq ~ I/wq = 0{1/N). This is shown in the right panel of Fig. 14. Let us thus define 
F'{z) = f' [z) + l/(z — Eo), where f' {z) contains the roots on the arc. Taking only one root away from 



22 



0.2 


-0.2 

0.2 


-0.2 

0.2 


-0.2 

0.2 


-0.2 





-0.2 





An^ ' An^-5 - 


— X 

- ,'*,X, X 




0.5 1 


-' k 'X 
- X 




- X 

-, ^ ,x 


m* • — • ■ ■ 


0.5 1 


- ' ^ x 

" X 




X 

-, ^ ,x 


• «t • — • ■ • 


0.5 1 


- ' X 'X 

- X 




- X 

- , X ,x 


• • «< » • • 


0.5 1 


r x' 


urn. • — • • • 


— X 

- X, ^ 





X 

X 




xr 




, 1 , 




- X 

X 




y 




-m — • 

, 1 , 











0.5 


1 


X 

-X 


X 




1 1 


II 




L X 


X 




1 1 


II 











5 


] 


p-i r- 

X 


-X- 




1 1 
• K « 




- 


_ X 










- 









0.5 


1 


- X ' 


-X 






1 1 1 1 




-X 

- ,>< , 


-X 




K» • 


• 

1 1 1 1 











0.5 


1 


- 1 X 
-X 




A 


1 1 


1 1 1 1 




-X 

- X 




4<- 


1 1 


-• 9^ 

1 1 , 1 





-0.2 

0.2 




0.2 


-0.2 

0.2 


-0.2 

0.2 


-0.2 



0.5 



0.5 



1 



FIG. 13: (color online) The location of those BA numbers which are 0(1) (red crosses) in the complex plane for 
the ground state and the lowest nine excitations in the sector with 5*^^ = (top left: ground state, bottom left: 
level 4) with Nb = 15, a;i = 2, _B = 2 in Eq. (2) for the couplings. The smallest six couplings are shown here (blue 
dots). Not shown is ujq = 0{Nb). 



the arc does not modify the arc in leading order. such that Fq{z) — Fo{z). This means that in leading 
order, the roots on the arc are decoupled from Eq- 

Let us now focus on the ground state of the central spin model with /ig = 0(1). If we assume that the 

roots Ek=i A/i, in the ground state of the central spin model are still described by an arc in the complex 

plane for large particle number, we have 71 = —F (0) + ^ + ^ — ho. This implies for the leading order 

(0) 

7i 

7i"^ ^ _Fo(0) + 1 + ^-/^0 (104) 



E 



2d 2 

yiT^A? 1 



^-V (105) 



where we have used the correspondence between BCS- and central spin parameters in Eqs. (11) and (47). 
Note that Aq — I/eq drops out in the first line. 

The two parameters 6, fi are now determined by the asymptotes of ^0(2:), 

ho = - lim Fo(z) = i y = (106) 



3= 



JVi, . 

ly — Ai 



2S^,, = - lim z [Foiz) + ho] = -_= + ^ , ' , ,,2 ' ^^^^'^ 

where we have set Eo ~ 0^ . The first of these equations coincides with Eq. (41). The second equation 
(107) is identical to Eq. (40). 
We have verified numerically that 7^°^ + ho = 0{l/Nb) for fields ho = 0(1). More generally, for 

^fot = C(l), the quantity 7^°^ + ho is of the order O{ho/Nb), as can be seen from a simple physical 
argument: Since the central spin is coupled to Nf, bath spins, for S^^^ = 0(1) it experiences an effective 
field /lo/iVb, and so does each bath spin. Thus the leading /iQ-dependent part of the spin-spin correlation 
function scales as ~ ho/N^, which yields a contribution ^ ho/Nf, to the energy. For S^^^. = 0(1) this 
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FIG. 14; a) Bethe roots ujk ~ 0{1), k = 1, . . . ,7 for < ho < 1.0. The field drives the roots away from the origin; 
the real root cjq, shown in b), moves towards the origin, c) The inverse numbers Ek = including fe = (on 

the real axis close to the origin, d) is a zoom showing Eo only). Here, the field makes the roots move towards the 
origin, except for Eq, which is shifted to larger real values. 



is just the leading contribution from 7^"'' + Hq. For a small central field, this has been demonstrated in 
Eqs. (80), (82) and (88). 

Since 71 is small compared to loq ~ 0{Nh), we can still use Eq. (69) to determine ujq iteratively. Thus 
Eq. (78) is still applicable for the energy, resulting now in 
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In the small-field limit, the results of the previous section are recovered. 

Comparing Eq. (109) with Eq. (42), one identifies the leading classical contribution due to from 

Eq. (42) stemming from the roots on the arc, i.e. "ff''^ . The root ojo encodes additional quantum-mechanical 
fluctuations which are of the same order of magnitude as the classical /ig-terms. 

In analogy to the energy, quantum fluctuations are also present in the correlation functions. In leading 
order, {Sj) is given by the classical expressions (38) and (39). Fluctuations are due to loq, which would 
yield a contribution ~ <P to (Sj), as in Eq. (67). However, the situation is different for the two-point 
function (Sq ■ Sj): The quantum fluctuations in the energy lead to contributions of order 0{d) in the 
two-point function, cf. Eq. (64). Taking together Eqs. (43) and (64), one obtains 



(5n • S 



2^1 



(0) 



dA, 



(110) 



Whereas Eq. (93) is valid in the weak-field regime ho < d only, Eq. (110) gives the field dependence and 
the leading finite-size effects also for stronger fields > d. This result is compared to numerical data 
from complete diagonalization in Fig. 5, showing very good agreement. Moreover, from Eq. (110), it is 
clear how to separate classical from quantum fluctuations, giving nice insight into the essential physics 
of the model. 
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V. CONCLUSION 



We have studied the exact solution of the central spin model, focussing on spectral properties and static 
correlators. In particular, it is possible to analyze the magnetization profile and the two-point correlation 
function using a classical approximation, exact diagonalization, and the Bethe ansatz solution as three 
independent methods. 

The exact magnetization profile of the quantum model follows the classical approximation very well 
already for small system sizes. For a given distribution of coupling parameters an increasing central 
field typically enhances the antiferromagnetic alignment of nearby spins, while it favours a ferromagnetic 
alignment with the outer spin. The total magnetization of the system is typically small. 

For the two-point correlation function a similar tendency can be observed, but the classical solution 
must be significantly corrected by quantum fiuctuation terms as given in Eq. (110). Only for the outermost 
spins the classical solution tends to become exact. The reason for this is that in all cases we considered, 
classical contributions arc encoded by the moments x^i of the couplings, whereas quantum fluctuations 
are expressed in terms of the moments xi, i > 0. This means that the outer region of the quantum 
dot, where the nuclear spins are coupled weakly to the electron spin, are governed by classical physics, 
whereas the inner region experiences stronger quantum fluctuations, due to the larger spin-exchange. 

The classical approach is analogous to the original BCS mean field solution of the superconducting 
state. Typically the classical approximation works better for the BCS model since quantum, i.e. finite- 
size contributions are sub-leading compared to the mean-field solution, whereas in the central spin model 
both can be of the same order in the central spin model, depending on the quantity under consideration. 
The reason for this is that the pairing amplitude g — 0{1/N) in the BCS model, whereas the analogous 
parameter Hq = g^^ = C'(l) in the central spin model. In view of tunable interactions in ultracold gases, 
this could lead to the possibility of a new pairing phase for attractive electrons with fixed particle number, 
when the attraction g is of order one. 

After having demonstrated how to obtain the classical contributions from the exact quantum- 
mechanical solution, we must emphasize that if Hq ^ 0, the expectation value {A\Sj'^\A) vanishes for 
all eigenstates A. This is necessarily so, since A must have a definite magnetization S^^^ unless there is 
an accidental degeneracy in the system. For the equivalent BCS model this means that the BCS order 
parameter (c]j^c^-f) is identically zero for finite quantum systems. Technically, the well-known sponta- 
neous symmetry breaking can therefore only be realized in the thermodynamic limit in the BCS model, 
despite the fact that a description in terms of the mean field solution (i.e. classical vectors rrij) gives 
quantitatively good results also for finite systems. This is in contrast to the prototypical example for 
symmetry breaking in ferromagnets, where the ground state and excited states generically already carry 
a non-zero expectation value of the order parameter for finite system sizes. 

Our results are of direct importance for the study of non-equilibrium dynamics: The understanding 
of the magnetization profile of eigenstates allows to estimate overlaps of eigenstates with those non- 
eigenstates which are realistic initial states in the time-evolution of the electron coupled to the nuclear 
spins. The computation of those overlaps is crucial in order to estimate the decohcrence time. We leave 
this as a promising route for future research here. 

More generally, the study of classical and quantum contributions during the time evolution of non- 
equilibrium dynamics remains an important open question for future research. 
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